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Local anomalies around the third peak in the CMB 
angular power spectrum of the WMAP 7- year data 
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ABSTRACT 

We estimate the power spectra of CMB temperature anisotropy in localized regions 
on the sky using the WMAP 7-year data. Here, we report that the north hat and 
the south hat regions at the high Galactic latitude ^ 30°) show anomaly in the 
power spectrum amplitude around the third peak, which is statistically significant up 
to 3a. We try to figure out the cause of the observed anomaly by analyzing the low 
Galactic latitude (\b\ < 30°) regions where the galaxy contamination is expected to be 
stronger, and regions that are weakly or strongly dominated by the WMAP instrument 
noise. We also consider the possible effect of unresolved radio point sources. We found 
another but less statistically significant anomaly in the low Galactic latitude north and 
south regions whose behaviour is opposite to the one at the high latitude. Our analysis 
shows that the observed north-south anomaly at high latitude becomes weaker on the 
regions with high number of observations (weak instrument noise) , suggesting that the 
anomaly is significant at sky regions that are dominated by the WMAP instrument 
noise. 
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1 INTRODUCTION 

The cosmic microwave background radiation (hereafter 
CMB) provides us with a wealth of information on the his- 
tory of the universe. The thermal black body nature of 
the CMB energy spectrum is now considered as the firm 
evidence of the hot big bang scenario for the beginning 
of the universe jAIpher fc Herman! Il948l ; iDicke et all 1 19651 ; 
iPenzias fc Wilson! 1 19651 ). The existence of large-scale struc- 
ture in the universe also implies that there were primordial 
density perturbations as the seeds for structure formation. 
It was expected that these inhomogeneities would have the 



and polarization w i th high resolution and sensitivity 
|Bennett et all 120031 ; Ijarosik et all 120111 ). For every data 
release, the WMAP team presented their estimation of 
the angular po wer spectra for temperature and polariza- 
tion anisotropy ( Hinshaw ct al. 2003, 2007: Page et alj|2003l . 
120071 : iNolta et all 120091 : lLarson et all l201lh . By compar- 



tions (anisotropy) ( 


Sachs , 


z Wolfdl 19671 


Bond & Efstathiou 


|1987|) 


The CMB 



~y3l97d : 



covered by the Cosmic Background Explorer (COBE) Dif- 
feren tial Microwave Radiometers experiment (|Smoot et al.l 
Il992l ) and has been confirmed b y many ground-based 
and balloon-borne e xperiments (see iHu fc Dodelsonl 120021 : 
IScott fc Smootll2006l for reviews and references therein). 

Recently, the Wilkinson Microwave Anisotropy Probe 
(WMAP) has opened a new window to the precision cos- 
mology by measuring the CMB temperature anisotropy 
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ing the measured CMB power spectra with the theoret- 
ical prediction, the WMAP team determine d the cosmo- 
logical parameters with a few % precis ion jSpergel et al.l 
I2003I . 120071 ; iKomatsu et~aH 120091 . l201ll ). and found that 
the observed CMB fluctuations are consistent with predic- 
tions of the concordance ACDM model with scale-invariant 
adiaba tic fluctuations generated during the inflationary 
epoch (|Spergel et al.ll2003l ; |Pe"iris et al.|[2003l : lKomatsu et all 
l201ll ). Recent ground-based and balloon-borne experiments 
that have performed the CMB power spectrum measure- 
ment and the cosmological p arameter estimation include the 
South Pole T elescope (SPT; iKeisler et alj|20lil ). the QUaD 
experiment |Brown et alj 12009] ) . Arcminute Cosmolog y 



Bolometer Array Receiver ( ACBAR; iReichardt et al.l 



the Cosmic Background Imager (CBI; iMason et al 



2008), 



2003), 



the Atacama Cosmology Telescope (ACT; iFowler et al.l 
20ld). the Degree A ngular Scale Inte rferometer (DAS I 



Carlstrom et al.ll2003l). BOOM ERANG (I Jones et al.N2006l) 
Arch eops jBenoit et all 120031 ). and MAXIMA (|Lee et all 
120011 ). 

In the CMB data analysis the two-point statistics such 
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as the correlation function and the power spectrum has 
been widely used. In particular, the relation between the 
CMB angular power spectrum and the cosmological physics 
is well understood, and the tight constraints on the cos- 
mological parameters can be directly obtained by compar- 
ing the measured power spectrum with the theoretical pre- 
diction. Therefore, the accurate estimation of the angular 
power spectrum from the observed CMB maps is the essen- 
tial step. The efficient techniques to measure the angular 
power spectrum from the CMB temperature fluctuations 



oped (e.g., Gorski 1994; iTegmarkl 1 19971; iBond et all 


1998; 


Oh et al. Il999l; ISzamidi et al. 


2001; 


Wandelt et al.l 


2001; 


Hansen et al. 2002; Hivon et al. 


20021; 


Mortlock et all 


2002; 


Hinshaw et al. 20031; Wandelt & Hansen 20031; Chon et al. 


20041; lEfstathioull2004l; lEriksen et alj|2004al; IWandelt et al. 


20041; iBrown et al.|2004|Polenta et al.ll2005; Fav et al.l 


200? ; 


Dahlen & Simons! 120081; iDas et al.ll2009l; iMitra et all 


200S; 


Ansari & Magnevillell2010l; IChiang & Chenll201lf). 



Until now, the analysis of the WMAP CMB data 
has been made using the whole sky area except for 



stronglv contaminated regions dHinshaw et al.l 


20031. 


2007|; 


Nolta et al.ll2009l; 


Larson et al .11 20 1 ll; 1 S aha et al. 


120061. 


2008; 


Souradeep et al. 


2006: Eriksen et al.l 2007a; Samal et alj 


20ld; iBasak & Delabrouilld l201lf). On the other hand. 



there was only a small number of studies on the 
power spectr um measurement on the partial regions of 



the sky (e.g., Eriksen et al. 


2004bl; Hansen et all l2004allbl; 


lAnsari & Magncvillc 20101; 


Chiang & Chenl 201 ll). In this 



work, we measure the angular power spectra from the 
WMAP 7-year temperature anisotropy data set on some 
specified regions of the sky. We found that at high Galactic 
latitude regions there is the north-south anomaly or asym- 
metry in the power amplitude around the third peak of 
the angular power spectrum, which is the main result of 
this paper. Our result differs from the well-known hemi- 
spherical asymmetry in the angular power sp ectrum and 
the genus topology at large angular scales JParkl |2004|; 
Eriksen et al.l I2004bll3. l2007bl; lHansen et al.l l2004bl ; iBernuil 
20081 ; lHansen et al.ll2009l ; iHoftuft et al.ll2009f ). 

This paper is organized as follows. In section 2, we de- 
scribe the method of how to measure the angular power 
spectrum. The application of the method to the WMAP 
7-year data and the comparison with the WMAP team's 
measurement are shown in section 3. In section 4, we mea- 
sure the power spectra on various sky regions defined with 
some criteria, and present our discovery of the north-south 
anomaly around the third peak of the angular power spec- 
trum at high Galactic latitude regions. We also discuss the 
potential origin of such a north-south anomaly. Section 5 is 
the conclusion of this work. T hroughout th i s work we h ave 
used the HEALPix software (|G6rski et all 1 19991 . 120051 ) to 
measure the pseudo angular power spectra and to generate 
the WMAP mock data sets. 



2 ANGULAR POWER SPECTRUM 
ESTIMATION METHOD 

In this section we briefly review how the angular power 
spectrum is measured from the observed CMB temperature 



anisotropy maps. Throughout this work we do not consider 
the CMB polarization. 

In the ideal situation where the temperature fluctua- 
tions on the whole sky are observed, we can expand the tem- 
perature distribution T(n) in terms of spherical harmonics 
as 



T(n) = £ £ a <™ F '™( n )> 



(1) 



1 — 2 m=-l 



where n denotes the angular position on the sky, Yj m (n) 
is the spherical harmonic basis function, and the monopole 
(I — 0) and dipole (I — 1) components have been neglected. 
From the orthogonality condition of the spherical harmonics, 
the coefficients a; m can be obtained from 



dfiT(n)F ; * m (n) 



C2) 



where dS} denotes the differential solid angle on the sky. If 
the temperature anisotropy is statistically isotropic, then the 
variance of the harmonic coefficients ai m is independent of m 
and the angular power spectrum C; is given as the ensemble 
average of two-point product of the harmonic coefficients, 



(ai m a Vm ,) — Ci<5;;'<5?7 



(3) 



where the bracket represents the ensemble average and 5 U i is 
the Kronecker delta symbol (8 u t — 1 for 1 = 1', and <5 ;; ' = 
for I 7^ I'). For a Gaussian temperature distribution, it is 
known that all the statistical information is included in the 
two-point statistics like the angular power spectrum C,. 

For the single realization of temperature fluctuations on 
the sky, the ensemble average is replaced with the simple av- 
erage of independent modes belonging to the same multipole 
I and the angular power spectrum becomes 



sky 



1 



21 + 1 



X 



"■lr, 



with uncertainty due to the cosmic variance 



AC, = 



21 + 1 



C, 



(4) 



(5) 



In the practical situation, however, we usually exclude 
some portion of the sky area where the contamination due 
to the Galactic emission and the strong radio point sources 
is expected to affect the statistics of the CMB anisotropy 
significantly. The limited angular resolution, the instrument 
noise, and the finite pixels of the processed map also pre- 
vent us from using the formulas, Eqs. |l}-([5|, which are 
valid only in the ideal situation. The incomplete sky cover- 
age with an arbitrary geometry and the limited instrumen- 
tal performance break the orthogonality relation between 
the spherical harmonic functions. Thus, measuring angular 
power spectrum from the realistic observational data is more 
complicated. 

There are two popular ways of measuring the angular 
power spectrum from the observed CMB data with incom- 
plete sky coverage. One is the maximum likelih ood estima- 
tion method based o n the Bayesian theorem (|Bond et alj 
1 19981 ; iTegmarkl 1 19971 ). and the other is the pseudo power 
spectrum method that applies the Fourie r or harmonic 
tran sformation d ir ectly (jHivon et al.l 120021 ). For the lat- 
ter, iHivon et all (|2002l ) developed the so called MAS- 
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TER (Monte Carlo Apodised Spherical Transform Esti- 
mator) method which estimates the angular power spec- 
trum from the high resolution CMB map data. For mega- 
pixelized data sets like those fr om WMAP or Planck 
jPlanck Collaboration et alj|2011ah . the MASTER method 
is faster and more efficient than the maximum likelihood 
method. 

We refer to the pseudo angular power spectrum as C; to 
distinguish from the true angular power spectrum C'i which 
we try to estimate. Two quantities differ from each other 
due to the incomplete sky coverage, the limited angular res- 
olution and pi xel size, and the instrument noise, an d they 
are related by (|Hansen et alj|2002l ; iHivon et al1l2002l l 



{Ci) = Y J M W F v Bl{C v ) + {N l ), 



(6) 



where Bi is the beam transfer function, Fi is the pixel trans- 
fer function due to finite pixel size of the CMB map, (Ni) 
is the average pseudo noise power spectrum due to the in- 
strument noise, and M ;; / is the mode-mode coupling matrix 
that incorporates all the effect due to the incomplete sky 
coverage. The mode-mode coupling matrix is written as 



Mi, i, = 



2l 2 + l 

47T 



h 









{3 




(7) 



where Wi is the power spectrum of the window function for 
incomplete sky coverage (see section [qT2"1i . and the last factor 
on the right-hand side enclosed with the parenthesis repre- 
sents the Wigner 3-j symbol. We use the Fortran language 
and the Mathematic jj software t o calcu late t he Wigner 3- 
j sym bols; see also IHivon et al.l (|2002T l and iBrown et al.l 
(|2005h for the numerical computation of the Wigner 3-j sym- 
bols. 

To reduce the statistical variance of the measured an- 
gular power spectrum due to the cosmic variance and the 
instrument noise, we usually a verage powers with in the ap- 
propriate I bins. According to IHivon et all (|2002l ). the bin- 
ning operator is defined as 



2tt ,(6+1) ,W ' 



if 9 < l (b ^ < 1 <r I 
11 z *s 'low ' 'low 

otherwise. 



(6+1) 



0, 



The reciprocal operator is defined as 

,(b) 



0, 



if 2 SS C < * 
otherwise, 



(b+l) 



where b is the bin index. For ease of comparison, we use 
exactly the same /-bi nning adopted by t he WMAP teams 
(b runs from 1 to 45; lLarson et al.ll201ll ). The true power 
spectrum can be estimated from 



C b = K-}P b n(C t - (N t )), 



where 



P bl M w F v Bt,( 



(8) 



(0) 



is a kernel matrix that describes the mode-mode couplings 
due to the survey geometry, the finite resolution and pixel 
size, and the binning process. 

Ideally the average pseudo noise power spectrum (Ni) 



can be estimated by the Monte Carlo simulation that mim- 
ics the instrument noise. However, because of the statistical 
fluctuations in the noise power spectra itself, it is rather dif- 
ficult to obtain an accurate estimation of the power contri- 
bution due to the instrument noise from the single observed 
data. In our analysis we evade this problem by measuring 
the cross-power spectra between different channel data sets, 
in the same way as adopted by the WMAP team (see the 
next section). 



3 APPLICATION TO THE WMAP 7- YEAR 
DATA 

3.1 WMAP 7-year data 

The WMAP mission was designed to make the full sky 
CMB temperature and polariza tion maps with high accu- 
racy, precision, and reliability |Bennett et al.l I2003T 1 . The 
WMAP instrument has 10 differencing assemblies (DAs) 
spanning five frequency bands from 23 to 94 GHz: one DA 
each at 23 GHz (Kl) and 33 GHz (Kal), two DAs each at 
41 GHz (Ql, Q2) and 61 GHz (VI, V2), and four DAs at 
94 GHz (Wl, W2, W3, W4), with 0?82, 0?62, 0?49, 0?33, 
and 0?21 FWHM beam widths, respectively. In this work 
we use the foreground-reduced WMAP 7-year temperature 
fluctuation maps that are prepared in the HEALPix format 
with Aside = 512 (resolution 9; r9). The total number of 
pixels of each map is 12 x A a 2 idc = 3, 145, 728. Following the 
WMAP team's analysis, we use only V and W band maps to 
reduce any possible contamination due to the Galactic fore- 
grounds. We use the WMAP beam transfer function (Bi) 
and the number of observations (A Q b s ) for each DA. All the 
data sets used in this work are available on the NASA's 
Legacy Archive for Microwave Background Data Analysis 
(LAMBDA0 



3.2 Weighting schemes 

The window function assigns a weight to the temperature 
fluctuation on each pixel before the harmonic transforma- 
tion is performed. In our analysis the window function maps 
have been produced based on two weighting schemes, the 
uniform weighting scheme an d the inverse-noise weighting 
scheme (|Hinshaw et al.l 120031 ). The WMAP team provides 
mask maps which assigns unity to pixels that are used in 
the analysis and zero to pixels that are excluded to avoid 
the foreground contamination. We use the KQ85 mask map 
with resolution 9 that includes about 78.3% of the whole 
sky area. The mask map excludes the Galactic plane region 
with the strong Galactic emission and the circular areas with 
0?6 radius centered on the strong radio point sources. The 
uniform weighting scheme uses the mask map M(p) as the 
window function directly, 



W(p) = M(p), 



(10) 



where M(p) denotes a mask map on a pixel p on the sky. 
The power spectrum at high I region is usually dominated 
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5000 10 4 1.5xl0 4 

Number of observations (W4 band) 

Figure 1. (Top) A map of number of observations N oos (p) at the 
WMAP 7- year W4 frequency channel (DA), multiplied with the 
KQ85 mask map M(p) [see Eq. II12I I1 , Dark blue color corresponds 
to zero values and represents regions that are excluded by the 
KQ85 mask map, while dark red color denotes the value exceeding 
N bs 6000. The Mollweide projection in Galactic coordinates 
is used to display this map, where the Galactic center is located 
at the center, the Galactic longitude increases from center to left 
(I = 0°-180°) and from right to center (I = 180°-360°), and 
the Galactic latitude increases from bottom to top (b = —90°- 
+90° ) . The regions with large number of observations correspond 
to the ecliptic pole regions. (Bottom) A histogram of number of 
observations at the same frequency channel. 



by the instrument noise. In the WMAP data the noise level 
on each pixel is modeled by 

(To 



a{p) = 



s/NabsJpj' 



(11) 



where iV bs (p) is the number of observations on the pixel p 
and ao is the global noise level (<to = 3.319, 2.955, 5.906, 
6.572, 6.941, 6.778 mK for VI, V2, Wl, W2, W3, W4 DAs, 
respectively). To reduce the effect of the instrument noise, 
the inverse-noise weighting scheme is defined as the product 
of the mask map and the map of number of observations, 



W(p) = M(p)N ohe (p). 



(12) 



The window function at 94 GHz W4 frequency channel (DA) 
is shown in Figure [1] together with a histogram of the num- 
ber of observations. Although the WMAP team applies the 
uniform weighting scheme for I < 600 and the inve rse-noise 
weighting scheme for I > 600 (|Larson et alj|20l"ll ). in this 
work we simply apply the single weighting scheme over the 
whole range of I considered (2 ^ Z ^ 1200) during the power 
spectrum estimation and compare the results based on the 
two weighting schemes. 

3.3 Angular power spectrum measured from the 
WMAP data 

In the WMAP team's data analysis, the power spectrum 
at low multipoles (I ^ 32) was measured by a Blackwell- 



Rao estimator that is applied to a chain of Gibbs samples 
obtained from the foreground-cleaned CMB map while the 
power spectrum at high multipoles (32 < I ^ 1200) by 
the MASTER pseu do power spectrum estimation method 
|Larson et al.l l201ll ). Throughout this work, the angular 
power spectrum over the whole I range (2 I ^ 1200) is 
measured based on the pseudo power spectrum estimation 
method because our primary attention is paid on the high 
I region where the peaks are located. Thus, as shown below 
our results at low multipoles are somewhat different from 
the WMAP team's results. 

The WMAP V and W bands have two (VI and V2) 
and four (Wl, W2, W3, and W4) separate frequency chan- 
nels (DAs). Using the Anafast program in the HEALPix 
package, we measured 15 pseudo cross power spectra (Ci) 
for different channel combinations (V1W1, V1W2, and so 
on) based on the uniform and the inverse-noise weighting 
schemes. During the production of each pseudo cross power 
spectrum, Anafast program estimates two separate sets of 
harmonic coefficients for each DA using the formula Eq. ([2]) 
where T(n) is now replaced with W (n)T(n) and W(n) is the 
window function defined in Eqs. (| 10p and (I12p for the corre- 
sponding DA. We set the maximum multipole as I — 1200. 

For each channel combination, we calculate the mode- 
mode coupling matrix M n i and the kernel matrix K bb i us- 
ing the beam transfer, the pixel transfer funct ions, and the 
same Z-binnhig as defined by the WMAP team (|Larson et al.l 
l201lfi . For the power spectrum of window function in Eq. (J7J) 
we use the cross power spectrum of window functions for the 
corresponding DA combination, which is obtained from the 
Anafast program by inserting the window functions at two 
different frequency channels as the input data. The squared 
beam transfer function B? appearing in Eq. ([6]) is also mod- 
ified into the product of beam transfer functions from the 
two different frequency channels. 

The final binned cross power spectrum Ct for each 
channel combination is obtained from Eq. @ neglecting 
the term for the average noise power spectrum (Ni). Since 
the properties of the instrument noise for each channel are 
independent, the process of cross-correlation statistically 
suppresses the contribution of the instrument noise to the 
power spectrum estimation. For each binned cross power 
spectrum we subtract the contribution due to the unre- 
solved radio point sources expected in the WMAP tempera- 
ture anisotropy m aps b y using the informa tion presented in 
iNoIta et all (20091 ) and [Larson et all |201ll ). The amplitude 
of the binned power spectrum due to the unresolved point 
sources can be written as 

Cps.b = ApsPbiS}. (13) 

Here Sf is the point-source spectral function given by 

V k V h i 



S\ = r{v h )r{v k t) 



Q 



P 



(14) 



where and vy denote the individual frequency channel 
belonging to i-th channel combination, r{v k ) is a conversion 
factor from antenna to thermodynamic temperature, i/q = 
40.7 GHz is the Q-band central f requency, 10 3 A VB = 9-0 ± 
0.7 fiK 2 , and /3 = -2 09 (see also iHuffenberger et al.1l2006l ; 
ISouradeep et al.ll2006l ). The subindex I is dummy because 
the spectral function does not depend on it. 
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To combine these cross power spectra into the single an- 
gular power spectrum, in fact we need the Monte Carlo sim- 
ulation data sets with the similar properties of the WMAP 
observational data. Using the Synfast program of HEALPix 
and assuming the concordance fiat ACDM model with pa- 
rameters, D. b h 2 = 0.02260, Q, c h 2 = 0.1123, fi A = 0.728, n s = 
0.963, t = 0. 087, Al(k = 002 M pc" 1 ) = 2.441 x 10~ 9 
(Table 14 of iKomatsu etafl l201lh . we have made 1000 
mock data sets that mimic the WMAP beam resolution 
and instrument noise for each channel. We use the CAMB 
software to obtain the t heoretical model power spectrum 
(ILewis fc Challinor1l2000l ). We assumed that the tempera- 
ture fluctuations and instrument noise follow the Gaussian 
distribution. In the simulation data sets we do not consider 
any effect due to the residual foreground contamination. We 
have analyzed the one thousand WMAP simulation mock 
data sets in the same way as the real data set is analyzed. 

Before combining the cross power spectra into the single 
power spectrum, each cross power spectrum has been aver- 
aged within /-bins specified to reduce the statistical fluctu- 
ations due to the cosmic variance and the instrument noise. 
We use t he same /-bins defi ned by the WMAP 7- year data 
analysis (|Larson et al.|[201ll ). Thus, in a case when we mea- 
sure the power spectrum from the WMAP data with a sky 
fraction much smaller than the area enclosed by KQ85 mask, 
the measured powers at low and high / regions become uncer- 
tain and fluctuate significantly with strong cross-correlations 
between adjacent bins (see below). We average the 15 binned 
cross power spectra into the single angular power spectrum 
by applying the combining algorithm which is similar to that 
used in the WMAP first year data analysis |Hinshaw et al.l 
2003). For each Z-bin, we construct a covariance matrix de- 
fined as 

(Efau)« = ([CI - Ct][Ci - C(}) + (E src )«, (15) 

where i and j denote the cross combination of the WMAP 
frequency channels (V1V2, V1W1, W1W3, and so on), b is 
the binning index (b = 1,. . . ,45), CI is the measured cross 
power spectrum for channel combination i, CI is the en- 
semble averaged cross power spectrum for the same channel 
combination expected in the concordance ACDM model (we 
have used one thousand WMAP mock data sets to estimate 
the ensemble averaged quantities), and 

(E src ); j = (PuSt)(P bl ,Sl,)aL (16) 

is the covariance matrix between errors expected to be 
caused during the subtraction of powers due to unresolved 
point sources. We use a 2 IC = (5A pB ) 2 = (0.0007 ^iK 2 sr) 2 . 
The final combined angular power spectrum is obtained by 

i=l < ~ , b\^{ u n)b , , 

Ub - ^15 ^15 / v -iyj ■ 

To estimate the uncertainty for the combined power spec- 
trum Cb, we obtain one thousand combined power spectra 
(Cj lm ) together with their average and standard deviation 
at each bin. The standard deviation measured at each bin is 
used as the uncertainty (or error bar) for the measured com- 
bined power spectrum. By comparing the theoretical model 
power spectrum with the average power spectrum obtained 
from the simulation data sets, we can check whether our 
power spectrum measurement algorithm works correctly or 
not. 



The left panel of Figure[2]shows the angular power spec- 
tra of the WMAP 7-year temperature maps measured with 
the uniform (red) and the inverse-noise (blue dots with error 
bars) weighting schemes. We have used the sky area defined 
by the KQ85 mask that was adopted in the WMAP team's 
temperature analysis. For a comparison, the WMAP team's 
measure ment has been show n together as grey dots with er- 
ror bars |Larson et al.|[201ll ) . Our estimations of the angular 
power spectrum in both the uniform and the inverse-noise 
weighting schemes are consistent with the WMAP team's es- 
timation. We note that the amplitude of the power spectrum 
at the 41th bin (851 ^ I sC 900) is slightly larger than the 
WMAP team's estimation, but they are statistically consis- 
tent with each other. Although statistically consistent, our 
result is a bit different from the WMAP team result at low 
I region because the WMAP team applied the different esti- 
mation method (Blackwell-Rao estimator based on Gibbs 
sampling) at this low / region while we simply used the 
pseudo power spectrum estimation method over the whole I 
range. 

The right panel of Figure [5] shows the averaged angu- 
lar power spectra measured from the 1000 mock data sets 
that mimics the WMAP instrument performance based on 
the concordance ACDM model. The averaged power spectra 
based on both the uniform and the inverse-noise weighting 
schemes coincide with each other and restore the assumed 
model power spectrum (green curve) within measurement 
uncertainties, which demonstrates that our measuring al- 
gorithm works correctly. However, our algorithm gives a 
slightly positive bias with respect to the true value at the 
last two bins at the highest multipoles, where the uncer- 
tainty due to the finite size of the pixels is significantly large. 
Comparing the results based on the two weighting schemes, 
we can see that in a case of the inverse-noise weight scheme 
the size of error bars is bigger (smaller) at lower (higher) 
multipoles. 



4 POWER SPECTRA OF VARIOUS LOCAL 
AREAS ON THE SKY 

Here we present the angular power spectra measured on 
various local regions on the sky defined based on the sim- 
ple criteria such as the Galactic latitude or the number-of- 
observations cuts. To define a partial sky region, we basically 
use the KQ85 mask map and the number of observations at 
the W band (W4 DA). By putting a limit on the Galactic 
latitude or the number of observations, we have obtained 
several local regions with different characters. 

Figure [3] summarizes the mask maps (or the window 
functions in the uniform weighting scheme) that are used in 
our power spectrum measurements. We estimate the angu- 
lar power spectra on each local area using the correspond- 
ing mask map as a window function. Especially, the north 
hat (b ^ 30°; b is the Galactic latitude) and the south hat 
(b ^ —30°) regions defined by the simple Galactic latitude 
cut show a difference in the power amplitude around the 
third peak of the angular power spectrum, which can be 
interpreted as the anomaly effect (see Fig. [4] below). To as- 
sess the statistical significance of the observed anomaly, we 
compare the difference of power spectrum amplitudes with 
the prediction of the fiducial flat ACDM model. Then, we 
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Figure 2. (Left) Angular power spectra measured from the WMAP 7-year temperature anisotropy maps. Our results are shown as red 
and blue curves with dots and error bars for uniform and inverse-noise weighting schemes, respectively. The sky area (78.27% of the 
whole sky) defined i n the KQ85 mask m ap has been used. For a comparison, the WMAP team's result is also shown as grey curve with 
dots and error bars iLarson et al.ll201 if) . (Right) Averaged a ngular power spectra measured from the one thousand WMAP mock data 
sets that are consistent with the concordance ACDM model l lKomatsu et al,ll201lh . Red color denotes the result based on the uniform 
weighting scheme while blue based on the inverse-noise weighting scheme. The power spectrum of the assumed concordance ACDM 
model is shown as green curve. 




(a) High latitude, north hat 
with b > 30° (24.11%) 




(e) North hat with high instru- 
ment noise (15.20%) 




(b) High latitude, south hat 
with 6 < -30° (23.98%) 




(f) South hat with high instru- 
ment noise (15.14%) 




(c) Low latitude, north region 
with 0° < b < 30° (15.15%) 




(g) North hat with low instru- 
ment noise (13.05%) 




(d) Low latitude, south reigion 
with -30° < b < 0° (15.03%) 




(h) South hat with low instru- 
ment noise (13.01%) 



Figure 3. Mask maps of various local areas defined by applying the Galactic latitude cuts [(a)-(d)] and the thresholds on the (smoothed) 
map of number-of-obscrvations [(e)— (h)] on the KQ85 mask map. The black area corresponds to M(p) = 1 while the grey area to M(p) = 0. 
Each number in the parenthesis indicates the fraction of sky area with M(p) = 1. 



search for the origin of the observed anomaly by considering 
possibilities due to the residual Galaxy contamination, the 
WMAP instrument noise, and unresolved point sources. 

4.1 North hat and south hat 

Figure[4]displays the angular power spectra measured on the 
north hat and the south hat regions. Each power spectrum 
has been estimated in both the uniform and the inverse- 
noise weighting schemes separately. The overall features in 
the measured angular power spectra are similar to the power 



spectrum measured on the whole sky region with KQ85 
mask (grey curves). However, the amplitudes of the power 
spectrum around the third peak show opposite deviations in 
the south hat and the north hat. 

We can see such a difference more dramatically in the 
middle panels of Figure [H where the difference of powers 
between the south hat (SH) and north hat (NH) are dis- 
played for an ease of comparison. We have also estimated 
the average and the variance of differences between SH and 
NH powers from the one thousand WMAP mock data sets, 
which are shown as dark green dots with la (dark grey) and 
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Figure 4. (Top) Angular power spectra measured on the north hat and the south hat regions based on the uniform (left) and the 
inverse- noise (right panel) weighting schemes. The power spectra measured on the north hat (b ^ 30°) are shown as red color while 
those on the south hat (6 —30°) as blue color. The grey color denotes the result measured on the whole sky enclosed by KQ85 mask 
map. (Middle) Difference of powers between the south hat (SH) and north hat (NH) are shown (orange color). The averaged difference 
estimated from the one thousand WMAP mock data sets are shown as dark green dots with lcr (dark grey) and 3cr (light grey) error 
bars. (Bottom) Histograms of differences between SH and NH powers in the 40th bin that corresponds to 801 ^ / ^ 850, measured from 
the ACDM based WMAP mock data sets, for uniform (left) and inverse-noise (right panel) weighting schemes. The vertical red dashed 
lines indicate the difference values measured from the WMAP 7-year data. 



3<r (light grey) error bars. Since we have assumed the homo- 
geneity and isotropy of the universe during the production of 
the WMAP mock data sets, the average of differences is ex- 
pected to be zero over the whole I range. We note that in the 
third peak (that corresponds to the 40th bin; 801 I ^ 850) 
the difference between SH and NH regions is statistically sig- 
nificant because it is located around the 3a deviating from 
the zero value. The histogram of differences between SH and 
NH powers in the 40th bin measured from the ACDM based 
WMAP mock data sets and the value measured from the 
WMAP data (vertical red dashed line) also confirms this 
result (bottom panels of Fig. [4| . 



4.2 Low latitude area 

The foreground contamination at high Galactic latitude re- 
gions like the north hat and the south hat is generally ex- 
pected to be small. Although the foreground model has been 
subtracted from the observed temperature fluctuations and 
the sky region that are significantly contaminated by the 
Galactic emission has been excluded by the KQ85 mask, we 
can expect that the residual foreground emission at low lati- 
tude regions may be stronger than high latitude regions such 
as the north hat and the south hat. We estimated angular 
power spectra on two separate regions at the low Galactic 
latitude defined as low latitude north (0° < b < 30°) and 
south (—30° < b < 0°) regions, where the Galactic plane 
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region has been excluded by the KQ85 mask. The results 
are presented in Figure [5] with a similar format as in Figure 

H 

The angular power spectra measured on the sky regions 
that are expected to be contaminated by the residual Galac- 
tic foreground emission show different features from those 
seen in the cases of the north and south hats. The behav- 
ior of power spectrum amplitudes around the third peak is 
opposite to the case of north versus south hat regions. The 
power on the low latitude north region is larger than that 
on the south region, and the difference between the north 
and the south regions is more significant at the 41th bin 
(851 ^ I ^ 900). The difference is again statistically sig- 
nificant with a 3<r deviation from the zero value (for uni- 
form weighting scheme), which is another anomaly found 
in this work. Although not shown here, if we consider the 
whole northern and southern hemispheres as two localized 
regions, then the anomalies at high and low Galactic lat- 
itudes are compensated and not noticeable. Furthermore, 
there are two more noticeable differences between the north 
and the south regions at the second bin (I — 3) and at the 
last bin (1101 ^ I ^ 1200). For the latter, the measured dif- 
ference exceeds 3a from the zero value for the inverse-noise 
weighting scheme (see also Table [T] below). Therefore, based 
on this comparison, it seems that the north-south anomaly 
around the third peak of the angular power spectrum ob- 
served on the high Galactic latitude north and south hat re- 
gions is not due to any possible residual Galactic foreground 
emission. 

4.3 Instrument noise 

Although the WMAP satellite probed the CMB temperature 
fluctuations on the whole sky, its scanning strategy is some- 
what inhomogeneous such that the regions near the ecliptic 
poles were probed many times as compared to the regions 
near the ecliptic plane. Thus, it is expected that the CMB 
anisotropy is strongly affected by the instrument noise in 
the ecliptic plane region. To quantify such an effect due to 
the WMAP instrument noise, we use the number of observa- 
tions, -/V bs(p), at the WMAP W4 frequency channel, which 
is shown in Figure [T] 

Defining regions with high and low instrument noise by 
simply putting threshold limits on the number of observa- 
tions results in the geometrical sky regions whose boundary 
is not smooth and the pixels near the boundary are not con- 
tiguous with many small islands outside the primary region, 
which prevents us to obtain the unbiased estimation of the 
power spectrum. To avoid this problem, we have used the 
Smoothing program of the HEALPix software to smooth the 
map of number-of-observations at the W4 frequency chan- 
nel with a Gaussian filter of 1° FWHM. After smoothing, 
we also exclude the island areas that are not included in the 
primary big regions. Then, we define the north hat (or the 
south hat) regions with high instrument noise by selecting 
pixels with N bs(p) < 2500 on the smoothed number-of- 
observations map. Similarly, to define the regions with low 
instrument noise, we set N \:, s (p) > 2035 on the same map. 
The threshold value for the number of observations at the 
region of high (low) instrument noise has been set to ob- 
tain the sky area of about 15% (13%) of the whole sky for 
-/Vobs(p) < 2500 (N bs{p) > 2035). For ecliptic plane regions 



p value (%) 

area Z-bin 







uniform 


inverse-noise 


High latitude (|6| > 30°) 


40 


0.1 


0.7 


— with high instrum. noise 


40 


0.5 


1.0 


— with low instrum. noise 


10 


31 


19 


Low latitude (|6| < 30°) 


41 


0.5 


5.0 




15 


8.5 


0.1 



Table 1. Probabilities of finding the north-south (absolute) dif- 
ference larger than the observed value, estimated based on one 
thousand ACDM based WMAP mock data sets. 

with high instrument noise (i.e., with low number of obser- 
vations), we choose the threshold value to include the larger 
sky area for the purpose of reducing the statistical variance 
in the power spectrum measurement. The sky regions de- 
fined in this way are shown in Figure 

We have measured angular power spectra on the ecliptic 
plane regions with small number of observations (defined as 
iVobs < 2500) that belong to the north hat or the south 
hat regions. The results are shown in Figure [6] The angular 
power spectra measured on the ecliptic plane regions (that 
belong to the north hat or the south hat regions) show the 
behavior that is similar to the case of the whole north hat 
and south hat regions (see Fig. [4]). The difference of power 
spectrum amplitudes between the south and the north hats 
still statistically significant, deviating from the zero value 
up to 3(7. 

We have also measured angular power spectra on the 
ecliptic pole regions with large number of observations (de- 
fined as N bs(p) > 2035), which are shown in Figure [7| 
Here, the power spectra measured on the ecliptic pole re- 
gions in the north hat and the south hat are consistent with 
each other in the I range around the third peak. They are 
also very consistent with the power spectrum measured on 
the whole sky with KQ85 mask (grey dots with error bars). 
Around the third peak, the power spectrum measured on 
the ecliptic pole regions in the north hat is very similar to 
the case of the whole north hat region, while the power spec- 
trum measured on the ecliptic pole regions in the south hat 
has decreased in amplitude compared with the case of the 
whole south hat region (see Fig. [4]). As a result, the observed 
north-south difference in the 40th bin of the power spec- 
trum is around la confidence limit and is not statistically 
significant any more (bottom panels of Fig. [7]). These results 
strongly suggest that the north-south anomaly around the 
third peak likely originates from the unknown systematic ef- 
fects contained on the sky regions that are affected by strong 
instrument noise. 

Based on the one thousand WMAP mock data sets, we 
count the number of cases where the magnitude of the north- 
south power difference is larger than the observed value for 
several chosen local regions. The results are summarized 
in Table [T] for particular /-bins where the high statistical 
significance is seen. The listed p-values are probabilities of 
finding the north-south (absolute) difference larger than the 
observed value. Although the statistical significance for the 
north-south anomaly at the high latitude regions becomes 
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Figure 5. Similar to Figure[4] but for regions at the low Galactic latitude areas (\b\ < 30°). The power spectra are measured on the 
north region (NR; 0° < b < 30°; red color) and on the south region (SR; —30° < b < 0°; blue color). In the bottom panels, shown are 
histograms of power differences at the 41th Z-bin (851 ^ I ^ 900) between south and north regions, expected in the concordance ACDM 
model, together with observed values from the WMAP 7-year data (vertical red dashed lines). 



weaker for the inverse-noise weighting scheme, the observed 
anomaly is still rare in the ACDM universe with p = 0.7%. 
Comparing the cases for high and low instrument noise at 
the high latitude regions strongly suggests that the north- 
south anomaly around the third peak come from the region 
dominated by the WMAP instrument noise. 

4.4 Point sources 

Recently, the Pla nck point source catalogue has been 
publicly available (|Planck Collaboration et al.l l2011bl ) . We 
take the Planck Early Release Compact Source Catalogue 
(ERCSC) at 100 GHz frequency band and exclude the 
sources at high Galactic latitude (|6| ^ 30°) whose angular 
position on the sky is located at the excluded region with 
M(p) — in the KQ85 mask map. The remaining 163 ra- 
dio point sources at |b| > 30° are considered as the sources 



that were newly discovered by the Planck satellite at 100 
GHz. By excluding (i.e., by assigning M(p) = at) the cir- 
cular areas centered on these point sources with radius 0?6 
in the WMAP team's KQ85 mask map, we have produced 
a new mask map (or window function) and measured the 
angular power spectra based on this new window function. 
Because the number of newly discovered sources are very 
small, the additional exclusion of point sources almost does 
not affect the angular power spectra on the north and the 
south hat regions. Therefore, at the present stage we cannot 
find any evidence that the north-south anomaly around the 
third peak of the angular power spectrum originates from 
the unresolved point sources. 
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Figure 6. Similar to Figure [4] but for regions with high instrument noise (iV D b s < 2500 on the number-of-obscrvations map at W4 DA 
smoothed with FWHM = 1°) that belong to the north hat (left) and the south hat (right). 



5 CONCLUSIONS 

In this work we have compared the angular power spectra 
measured on various local sky regions using the WMAP 7- 
year temperature anisotropy data. The sky regions stud- 
ied in this work are the low Galactic contamination regions 
at high Galactic latitude (north hat and south hat), the 
strong Galactic contamination regions at low Galactic lati- 
tude north and south regions, the regions dominated by the 
WMAP instrument noise (ecliptic plane regions), and the 
regions of low instrument noise (ecliptic pole regions) . 

We found that the power spectra around the third 
peak measured on the north hat and the south hat regions 
show an anomaly which is statistically significant, deviating 
around 3a from the ACDM model prediction. We have tried 
to identify the cause of this anomaly by performing the sim- 
ilar analysis on the low latitude regions and the regions with 
high or low instrument noise. Curiously, in the low Galactic 
latitude (|6| < 30°) there appears another but less statisti- 



cally significant north-south anomaly around the third peak, 
whose behaviour is opposite to the one seen in the high lat- 
itude regions and compensates the anomalies in the whole 
northern and southern hemispheres. At the present situa- 
tion, we cannot draw any firm conclusion for the origin of 
the observed anomaly. However, we found that the observed 
north-south anomaly maintains with the high statistical sig- 
nificance in the power spectra measured on the regions with 
high instrument noise, and the anomaly becomes weaker in 
the power spectra on the regions with low instrument noise. 
Thus, in our present analysis the observed anomaly is signif- 
icant on the sky regions that are dominated by the WMAP 
instrument noise. 

The location of the third peak in the angular power 
spectrum corresponds to the angular scale that approaches 
the WMAP resolution limit and that is dominated by 
the WMAP instrument noise. Because the Planck satel- 
lite probes the CMB anisotropy with higher angular res- 
olution and sensitivity than the WMAP, we expect that 
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Figure 7. Similar to Figure [4] but for regions with the low instrument noise (AT t, s > 2035) belonging to the north hat (left) and the 
south hat (right). 



the origin of the observed anomaly will be identified in 
more detail when the Planck data b ecomes publicly available 
(|Planck Collaboration et alj|2011ah . 
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